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Abstract 



We study the local semicircle law for Gaussian /3-ensembles at the edge 
of the spectrum. We prove that at the almost optimal level of rC'^^^'^'^ , the 
local semicircle law holds for all /? > 1 at the edge. The proof of the main 
theorem relies on the calculation of the moments of the tridiagonal model 
of Gaussian /^-ensembles up to the p„-moment where p„ = 0{n^^^~'). 
The result is the analogous to the result of Sinai and Soshnikov [U for 
Wigner matrices, but the combinatorics involved in the calculations are 
different. 



1 Introduction and Summary 



Our goal in this paper is to prove a local semicircle law at the level of 0{n^^"^^'^) 
near the spectral edges of the Gaussian /3-ensembles. The corresponding prob- 
lem for Wigner matrices at the edge has been established by Sinai and Soshnikov 
[Ti] and in the bulk by Erdos, Schlein and Yau 0; [H]- 

Definition 1. Let /3 > 1, the Gaussian ^-ensemble is an ensemble in R that 
have the following probability density: 

d/(Ai,...,A„) = G„,/3 n |A,-Ajf exp(-^^A?)I]dA, (1) 

l<2<j<n 1 i 

where Gn,i3 is a normalization constant. 
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For j3 = 1,2,4, this probability distribution corresponds to the eigenvalues of 
GOE, GUE and GSE respectively. Other known and related results regard- 
ing the /3-ensembles will be discussed in the next section. Using the moments 
method, one can establish the Wigner semicircle law for the /3-ensembles, e.g. 
H: 

Suppose Xi,i — 1, . . . , n follow the distribution of a /3-ensemble. Consider the 
rescaled eigenvalues Ai — and the empirical distribution function: 

Nn{X) = : Afe < A} 

n 

We have 

.A 

lim iV„(A) / p(u)du 



p(u) = 



where 

for u > 1 or u < — 1, 

-Vl — u-^ for — 1 < u < 1. 



Take an r„-neighborhood On of right end of the spectrum A = 1, where r„ = 
gQjj^g e > 0. Normalizing A^ = 1 — 6'fer„ and 

nrn 

at each 9k, we obtain a measure /i„ on the real line such that /in(M) — 

The main result of this paper is the proof of the following local semicircle law 

at the edge of the spectrum: 

Theorem 1. As n — >■ oo, /i„ converge weakly in probability on each finite 
interval to a measure fi concentrated on and almost continuous with respect 
to the Lesbesgue measure: 

f ifx>0 
dx [0 otherwise 

If fn > n'^~^^^ for some e > 0, then the measures converge vaguely to p with 
probability 1. 

Remark 1. The reason we called the above theorem a local semicircle law is 
that under the semicircle law, one can show that for a neighborhood of 1 of 
size of some fixed number S, the proportion of normalized eigenvalues in that 
neighborhood converges to the corresponding integral of the semicircle law. The 
main theorem states that the size of this neighborhood can be much finer, at 
the scale of r^^2/3+e convergence will still hold. In other words, for a 

neigborhood of size r„ near 1, there are ilnrf/"^ eigenvalues. If one rescales the 
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total number in this interval by nrn , this number with converge in probability 
to the appropriate integral of the semicircle law. Note that one would not expect 
the semicircle shape to hold at the edge when rn — n^"^^^ . In this sense, the 
local semicircle law above is almost optimal. 



The theorem above wih be proven using the moment method, similar to the 
method by Sinai and Soshnikov |14) : 

Theorem 2. Let An be a random matrix whose eigenvalues are given by equa- 
tion (Qp scaled by a factor of l/2y/nf3. Let pn ^ oo as n ~> od and pn = 
Q^^2/3-e^ /or some e > 0, then 

IE( jy^Pn) ^ I 23/2n(7rp3) ^^^(1 + 0(1)) if Pn is even, 
" [0, if Pn is odd. 



However we shaU see in the proof that the combinatorics involved is different 
and simpler in our case. The derivation of theorem [T] from theorem[5]is the same 
as in |14| and will be given in the appendix. It should also be remarked that 
the result in this paper is slightly weaker than that in [14 , as in their paper, 
only r„ << n'^^^ is required. 

The structure of the paper is the following. We first introduce a matrix model, 
proposed by Dumitriu and Edelman [5], for the /3-ensembles and recall some 
related results. We shall then prove theorem [5] in section [31 

For the rest of the paper, we shall use the following notion of an event depending 
on some index n having overwhelming probability: 

Definition 2. We say an event E holds with overwhelming probability if for all 
n, P(£') > 1 — Oc{n^'^) for every constant C . 



It should be observed that a union of events of overwhelming probability for 
some fixed k still holds with overwhelming probability. 



2 Gaussian Beta Ensembles 



As mentioned above, the special values of /3 = 1,2,4 corresponds to eigenval- 
ues of GOE, GUE and GSE respectively and there is an extensive literature on 
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them. For the edge of the spectrum, Tracy and Widom [15] , [12] has been able 
to estabhsh that, upon rescahng and centering, the top eigenvalue distribution 
converges to what is known as the Tracy- Widom distribution. Later Soshnikov 
[13] proved the edge universality for Wigner matrices. However, the case of 
general (3 remains unknown until the recent paper by Ramirez, Rider and Virag 
[T5] . One of the obstacles for studying /3-ensembles previously is the lack of 
matrix formulation for them. In their paper [5], Dumitriu and Edelman suc- 
ceeded in writing down a tridiagonal random matrix model whose eigenvalue 
distribution is given by equation ([1]): 



Theorem 3 (Dumitriu-Edelman) . Consider the matrix given by 



X(n-l)^ 





A: 



X{n-l)fi 
iV(0,2) 

X(n-2)I3 





7V(0,2) 





X(n-3)/3 



... X3^ N{Q,2) X2P 

X2P N{0,2) x& 

V xp ^(0,2) j 

(2) 

where iV(0, 2) denotes a random variable whose distribution follows the Gaussian 
distribution with mean and variance 2; Xk is a random variable having a chi 
distribution with k degree of freedom. The upper triangular part of the matrix 
consists of independent random variables and the matrix is symmetric. The 
eigenvalues of this matrix follows the distribution of the f3-ensemble. 



For GOE, the theorem can be easily proven by successive Householder trans- 
formation and this was observed by Trotter 17 . For general (3, the idea of the 
proof is to compute the Jacobian of the diagonization of the matrix of the form 
©. 

Another important result is the one by Ramirez, Rider and Virag |13| . whose 
work is based on previous work by Edelman and Sutton [6^. In their pa- 
per, Ramirez, Rider and Virag proved that the edge scaling parameter is n'^^^ 
and computed the asymptotics of the rescaled limiting distribution of the j3- 
ensembles: 

Theorem 4 (Ramirez, Rider, Virag). Let Xi, . . . , Xk be the k largest eigenvalue 
of the /3-ensemble, in descending order, and 

2 



where b' denote a white noise. Let Ai, . . . , be the k lowest eigenvalues ofHp, 
then 




2nl3) 



4 



converges in distribution to (Ai,...,Afc) Moreover the asymptotics of Ai are 
given by: 

P(Ai >a)= exp(-^/3a3/2(i + ^(i))) 

and 

P(Ai < -a) = exp(-^/3a3(l + o(l))) 



The proof of the above theorem is based on the observation that the tridiagonal 
matrix model can be viewed as a discretization of certain stochastic differential 
operator Hp and the probability distribution of the largest eigenvalue of the 
tridiagonal matrix will converge to the probability distribution of the eigenvalues 
associated stochastic differential operator as n — >■ oo. This point process is given 
by the so-called Airyp process and upon rescaling will converge to /i in our 
main theorem. Notice that if we consider the operator "Hoc, which is the free 
Laplacian, then the spectral measure of the operator under proper rescaling is 
H. 

Of course, these two results are only the tip of the iceberg of the vast field of 
random matrices, there are excellent monographs , e.g. [T], [TUl that discuss the 
latest developments in details. 

We see from the above theorem that, at the scaling of n^^'^, the limiting dis- 
tributions are different for each /3 at the edge of the spectrum. On the other 
hand, the Wigner semicircle law is universal for all /?, one can therefore ask 
the question at which scale does the edge look different for different /3, which 
leads us to consider the problem at hand. Moreover, one can also ask, given the 
availability of a matrix representation, whether the classical moment method 
will be able to give us information at the edge of the spectrum. It is also part 
of the goal of this paper to explore this possibility. Lastly, a local semicircle law 
is always of interest. The power ofa local semicircle law can be seen in recent 
papers of Erdos, Schlein and Yau 0, [S] and in Tao and Vu [16] in their proof 
of universality in Wigner ensembles and also in the recent paper of Bourgarde, 
Erdos and Yau [3]. Bao and Su [5] has also established through the Brownian 
Carousel a local semicircle law in the bulk for the Gaussian /3-ensemble. Also, 
one does not expect to prove a local result at the bulk using the moments as 
higher moments highlight the eigenvalues at the edge. 
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3 Combinatorial considerations 



Our goal is to understand the trace of higher degrees of the matrix An^p. For 
simphcity of notations, we consider the matrix scaled by I/a/Sj so that the 
off diagonal entries have variance n,n — 1, . . . , 1 and the diagonal entries have 
variance 2//3. The only difference is that we have to normalize the resulting 
calculations of moments by 2P"nP^/^. The traces are given by 



r 

where V is the set of indices {io, ik-i} 

It is customary to encode this information as a graph and to study the combi- 
natorics associated with this graph: 




12 3 n-2 n-1 n 



Figure 1 : Graph representation of a tridiagonal matrix 

We first note an easy simplification of our problem at hand, due to Dumitriu 
(see [3] for the proof): 

Lemma 1. E(Tryl,';) = nE(A^)ii 

The plan now is to first calculate the contributions from paths that do not go 
through self loops, i.e. paths where ij ^ ij+i for all j; after that we shall prove 
that the other paths contribute asymptotically lower order terms to the sum. 

Let us discuss of the odd moments first. Even though the offdiagonal entries 
are not symmetric distributions and have nonzero odd moments. One notices 
that the coefficients of the characteristics equation for the matrix involves only 
even powers of the offdiagonal entries, so we can replace the x-distribution by 
a symmetric distribution that has the same even moments as the x-distribution 
and the joint distribution of eigenvalues of the resulting matrix will be the same 
as that of the original. Since the distribution of the trace depends only on that 
of the eigenvalues, we see that the odd moments vanish. If f{x) is the density 
of the x-distribution, then the distribution g{x) equals to /(x)/2 for a; > and 
/(— x)/2 for a; < will be the desirable replacement. 
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We now return to the calculations of the even moments. Associate with each 
path on the graph a Dyck path of length 2k, where travelling to the right 
corresponds to an up edge and travelling to the left corresponds to a down 
edge. Throughout the rest of the paper, we shall most of the time think of Dyck 
paths as probabilistics objects, i.e. samples from a random walk X of length 2k 
conditioned on Xq = X2k = and Xm > for < m < 2fc. 

Before we start the proof in earnest, let us first take a step back and recall some 
relevent known facts and see what could cause problems: 

1) The number of Dyck paths of length 2k is given by the fc-th Catalan number 
Cfc. 2) The 2l-th moment of a chi-distributed random variable of degree with 
freedom k is given by rii=o(''' ~^ ■^*)- 

Each path in the sum ^ will contribute JliLo("' + J«) ^^r some ji. The hope 
is that asymptotically, for k — o{in?^^), the product will be n^{l + o(l)) and 
the worry is that some of the ji gets too small (negative) or too big. This 
corresponds to the case where the path goes too far to the right or stay on the 
same edge on the graph for too many steps. Therfore our goal is to show that 
there are very few of these problematic paths that they do not affect the total 
sum as n goes to infinity. 

The first step is to connect Dyck paths with a true random walk, so that we 
can use many of the known results and nice properties of the latter. We start 
with the following well-known lemma, which establishes the relationship between 
Dyck paths and Bernoulli bridges. 

Lemma 2. There is a one-to-one correspondance between + 1} x 

{Dyck paths of length 2k} -H- {Bernoulli Bridges of length 2k} 



Proof. The proof is to provide an explicit map. Attach a down path to the end 
of a Dyck path, resulting in a total of /c + 1 down paths and k up paths. Pick one 
of the fc + 1 down paths and make a cut before the down path. The Bernoulli 
bridge (from (0,0) to {2k, —1) with first path always downwards) is contructed 
by moving the second half to the origin and attaching the first half to it; 

□ 



The next lemma, relating the probability of an event of a Bernoulli bridge with 
one of a random walk, is first proven by Khorunzhiy and Marckert llj: 
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Figure 2: A sample Dyck path 




Figure 3: The BernoulU bridge corresponding to the above Dyck path cutting 
at the second down edge 

Lemma 3 (Khorunzhiy, Marckert). Let Sk = Xi + . . . + Xk be a simple random 

walk. Let Ti be the measurable events generated by Xi , . . . , X; and let P]" be the 
probability measures induced on the paths by Xi, . . . ,Xi. Let Fi be the measure 
Ff conditioned on Xq = and Si = —1. Then for any Tk- measureable event 

P^,+i(^fe)<CP^fe+i(Afe) = CP^(Afe) (4) 
for some constant C independent of n. 



Proof. The second equality is obvious as is J^/s-measurable. Suppose A is 
J^/j-measurable, then 

V2k+i{A\Sk = m,S2k+i = -1) = V2k+i{A\Sk = m) = V^{A\Sk = m), 

thus, 

Fl,^,{A)=F^,^,{A\S2k+i = -l) 

= Y,Ft{A\Sk = m)F^k+iiSk = m\S2k+i = -!)• 
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On the other hand, 



(k\(k+l^ 

Vm/ Vrn+1/ 
(2k+l\ 
Vm+i; 



ok f2k+l\ 



where Co := supj, sup„ (\k+r\ is finite. It follows that 



.+1 {A) < Co ^ PJ^ (A| 5fe = m)Pfe {Sk=Tn) 

m 



□ 



As discussed earlier, we would like to show that the probability of crossing too 
far to the right of the graph is small and this is indeed the case: 

Proposition 1. Let S2k be a random walk conditioned on So — S2k = and 

Si>0 for allO<l< 2k. Then for all k 

P( max Sk > < Cexp(-cX^) (5) 

0<l<2k 

for some constants c and C . 



Proof. Using the notations from the previous lemma, denote the Bernoulli ex- 
cursion of length 2k, the Bernoulli bridge from (0, 0) to (2fc, 0), and the random 
walk of length 2k by 82^, 5*2^., 5^ respectively. 

By lemma [21 we have 

max Sf — min Sf — 1 < max sf < max Sf — min Sf + 1 

0<l<2k 0<l<2k 0<l<2k 0</<2fc 0</<2fc 

So it suffices to show that 

P( max St - min St > Xk^''^) < Cexpf-cA^) 

Moreover, if 

max S\ — min S\ > \k^^^, 

0<l<2k 0<l<2k 
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then either 



max St - min St > ^Afc^/^ (6) 
o<;<fc o<i<k 2 



or 

max st - min 5?. > -Xk^/^. 

k+l<l<2k k+l<l<2k 2 

By reversing the Bernoulli bridge, it suffices to show the first inequality. The 
event satisfying the inequality ^ is contained in the union of the events 

max St > and min St < — fc^/^ 

0</<2fe 4 0<i<2fc 4 



By lemma [3l 



( max Si > < CoP( max S^ > -^fc^/^) 



< 2CoP(5fe > (7) 

where the second inequality is implied by the reflection principle. Hoeffding's 
inequality then implies the desired inequality: 

P( max St > — < CexY>(-cX^) 

□ 



Next, we shall consider the potential problem of a path staying on the same 
edge too many times. Again, the probability of this event is very small: 

Proposition 2. Let Si,0 < I < 2k be the Bernoulli excursion as above. Let Ti 
denote the number of I such that Si = i, for i — 1, . . . , 2k, then 

P(maxT, > fci/2+^) < Cexp(-cfc2^') 

i 

for any e > e' > and some constants C, c independent of k, X. 



Proof. As in the previous proposition, we shall connect the quantity in question 
with a similar quantity of a random walk. Clearly, since by mapping Dyck 
paths to Bernoulli bridges, the times one visits a level is mapped to at most two 
different levels, we have 

P(maxi;'= > < P(2max7;'' > Xk^/^) 

i i 

and as before we split the bridge into two halves 

P(2 maxT^'' > Xk^^^) < P( max Tj' > -k^'^) + P( max T^" > -k^'^) 

^ j * ~ ' ~ ^Q<i<k ' ~ 4 ^k+l<i<2k * ~ 4 
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By lemma [31 it suffices to show that 

P( max T,"" > < Ccxp{~cX^) 

0<i<k 4 

We wiU in fact first show that, for each —k<0<k, 

PiTr > ^fc'^') < Cexp(-cA2) 

Let q2k,r be the probability that a random walk of length 2k returns to exactly 
r times. It is well known (e.g. |T2]) that 

_ 1 /2fc-r 

q2k,r - 22/c-r I ^ 

We can then approximate the probability that there are less than returns 
to for a random walk of length 2k using Stirling's formula 



[2 

"(there are less than fc^/^^^ returns to 0) ~ -i/— / 

V 71" Jo 



J,l/2 + e 



For levels other than 0, we can condition on the first hitting time and similarly 
show that 

¥{T, > k^/^+') < Cex.j>{-ck^') 
Therefore, we have the bound 

P(maxT, > fc^/2+") < Cexp(-cfc2'') 

i 

for any e' < e □ 



With the two propositions above, we are ready to compute the contribution by 
the main paths. First of all the number of Dyck paths of length 2k is given 
by the Catalan number Ck = khi^k)- ^^^^ shown that jr^(^fc)(l - o(l)) 
of these paths have 'good properties', i.e. they are not in the exceptional set 
in either propositions. The contribution of each paths is of order n'^(l + o(l)) 
when k = o(rt^/"^). To see this, we simply compute the contribution 

(n + ii) .. . [n + ik) 

for that particular path. Not being in the exceptional set in the first proposition 
means that ij > —n}/^~^ for any (5 > and not being in the exceptional set 
in the second proposition means that ij < n}/^~^ for any 5 > Q. Therefore 
^3 — and the claim follows. 

It remains to show that the other paths do not contribute asymptotically to 
the sum. Let us first tackle the Dyck paths that do not have 'good properties'. 
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Consider the paths where maxj Tj is of the order g„ where qn = pl/'^~^^ for some 
I > 0, each of the paths will contriubute 0{n'P"{l + exp(n'))). By the previous 
proposition, the probability density of these paths are asymptotically of the 
order 0(exp(— Integrating over these paths, we can conclude that they 
are asymptotically of 0(71^"). 

Lastly we compute the paths where there are at least one loop, i.e., on the path 
{1, ii, 12, • . • , ipt^i, 1}, there exists / such that ii = ii+i: 

Proposition 3. The contribution to the sum from paths with 2k loops, i.e. 
paths {l,ii,i2, ... ,ip„~i,l} where there exists h, ... ,l2k such that ii. — ii.^i 
}orj = l,...,2kiso}o{Cp^/2n'^'^') 

Proof. First notice that the number of loops must be even as the odds moments 
of a Gaussian random variable are zero. As in proposition ([2]), we let Tj denote 
the number of times the path stays on level i. The number of paths with 2k 
loops are bounded by 




since the loops have to be paired on the 'same level' and removing the loops 
give a Dyck path of length j3„ — 2k. With overwhelming probability, each path 
contributes 0{n^~^ Y\-{2ai — l)\\). Indeed, if we remove the 2k loops, it reduces 
to a Dyck path that the earlier analysis applies, hence the factor n^~^ and each 
loop on the same level contributes the 2ai-th moment of the Gaussian random 
variable, which is of the order {2ai — 1)!!. For a given fc, the contribution 
from the terms in equation ([5]) is dominated by the term where I — k and 
ai — . . . — ai = 2 (one way to see this is that this is a similar argument as the 
one which shows that trace of moments of Wigner matrices are dominated by 
paths that traverse each edge twice after considering only the contribution from 
the loops, see e.g. [M], the key here being the moments of Gaussian does not 
grow too fast; another way is to calculate the contribution when I < k. This is 
bounded by (" ^ )/'^~'n^'^/'^~'^'^([j] !!)', the binomial term coming from picking 
I levels for the loops, the second term being each level must be filled by at least 
two loops, and the remaining loops can be anywhere, the third term is just 
the maximum possible of (^')(2ai — 1)!! and as explained below, maxi Ti < 

Summing over all Z < fc yields the desired dominance.) and is bounded 
by fik/a-ke^ik/s-ke _ 0(n^/3-«)^ since we have max^ Tj = 0{n^/^~'^) with 

overwhelming probability and therefore (^') = 0(rt'^'/'^~'^) with overwhelming 
probability and moreover, by proposition ([1}, the number of different to 
choose from in the sum is of 0{n^^^~'^) with overwhelming probability and so 
the number of choice is bounded by n^/^-'"^ . The exceptional paths can be dealt 
with similar to the ones without loops. Using the trivial bound of Cp^/2-k < C'p„ 
result in the proposition. □ 
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If we sum over all fc > 2 in the proposition above, we found that the contribution 
from the paths with loops is of the order 0(Cp,^/2nP"/^^'^). 

Combining all of the above, we have that E(A^^)ii is given by 

when pn is even and when pn is odd. Here Ck denotes the k-th Catalan 
number, which counts the number of Dyck paths between and 2k. It is well 
known that: 

and the result follows from lemma [T] after we rescaled by 2^'"nP"/^. 



4 Conclusion 



In this article we prove a local semicircle law at the edge of the spectrum for 
/3-ensembles. The proof relies on calculations of moments and can be applied 
to other tridiagonal matrices satisfying fairly mild conditions: that the diagonal 
entries remain 0(1) and the growth of the moments of the off diagonal entries 
is not too fast. It would be interesting to consider the same question for the 
^-Laguerre ensembles. 
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6 Appendix 



We shall provide the derivation of theorem [T] from theorem [2] for the readers' 
convenience. 
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Proof. Proof of theorem [T] To prove the theorem, it suffices to estabhsh the 
convergence of the Laplace transform of the measures /i„ to that of /i (see e.g. [9] 
for the relationship between convergence of Laplace transform and convergence 
in distribution): 

Let Sn — cr~^/2 and p„ = 2s„. We claim that the sequence n~^r^^^TrA^''" 
converges in probability to 

hm 4^E(TrA^-0 - Hm " - ^ 

This follows from Chebyshev's inequality once we show that the variance of 
n'^rn^^TrA^'^'^ is of o(l). This is true as: 

var(n-VJ^/2Tr^2s„-) ^ n-^r-^vaTiTiA^"") 

By theorem [21 the last quantity is of o(l). Similarly, n~^r„ '^^^TrA^''"+^ con- 
verges in probability to zero. Let 

A^. = 1 - 0krnXk > (9) 

Aj = -1 + T,r„Aj < (10) 

Then by definition, we have 

TrA^^" = ^(1 - r„0fc)'»^/'^'-"' + r.r^ f^^^/'^''"'^ (11) 

k j 

and 

TrA2-"+i = ^(1 - r„0fc)2((-/2)-.-')+i - ^(1 - r,r„)2«=/2)'-"")+i (12) 

k j 

If l^**:! < T'n ^^'^ and ItjI < rn^^^, the corresponding terms in pT|) and (IT2t are 
e~°^*(l + o(ry'^)) and e^^'^^ (1 + o(ry'^)). On the other hand, using the estimate 
of the expectation of TrA***" , the subsums in (ITT|) and ((T^) over 6k and Tj such 
that \9k\ > rl/^ and \Tj\ > r^'^ converge to zero in probability. The same holds 
for the subsums over these 9k and tj of J^k and J2j e"'^^^ . Therefore, we 

have 

(TrA^^" + TrA^'^'^+i - 2 ^ 6^'=^'=)^— -> 



3/2 
k 



and 



3/2 V 7rr3 
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If Tn = n ^ where 7 < 2/3, then it foUows from theorem [2] that the normaUzed 
traces n^^r„ ^^^TrA^''" and n~^7'„'^^^Tryl^''"+^ converges to nonrandom Umits 
with probabihty 1 , and hence the convergence of the measures /Xn in the main 
theorem also occurs with probabhty 1. □ 
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